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The physical processes that control the partition of released magnetic en¬ 
ergy between electrons and ions during reconnection is explored through particle- 
in-cell simulations and analytical techniques. We demonstrate that the de¬ 
velopment of a large-scale parallel electric held and its associated potential 
controls the relative heating of electrons and ions. The potential develops 
to restrain heated exhaust electrons and enhances their heating by conhn- 
ing electrons in the region where magnetic energy is released. Simultaneously 
the potential slows ions entering the exhaust below the Alfvenic speed ex¬ 
pected from the traditional counterstreaming picture of ion heating. Unex¬ 
pectedly, the magnitude of the potential and therefore the relative partition 
of energy between electrons and ions is not a constant but rather depends 
on the upstream parameters and specihcally the upstream electron normal¬ 
ized temperature (electron beta). These hndings suggest that the fraction 
of magnetic energy converted into the total thermal energy may be indepen¬ 
dent of upstream parameters. 
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1. Introduction 

Magnetic reconnection is a nniversal plasma process which converts stored magnetic 
energy into particle energy. The process is believed to be important in many astrophysical, 
solar, geophysical, and laboratory contexts. A principle topic in reconnection physics is 
the mechanism by which magnetic energy is partitioned into electron and ion thermal 
energy. A measure of this partition is the relative fraction of the available magnetic 
energy per particle W = (47r mi Uup) = m^ (or its asymmetric generalization 

mic\ asym [Phau ct tti, 2013; Shay et ai, 2014]) that goes to each class of particle; the 
supscript ”up” denotes the upstream value and Bmp is the reconnecting component of the 
magnetic held. 

Ion thermal energy often makes up a large fraction of the released magnetic energy 
during magnetic reconnection in both in the magnetosphere [Eastwood et ai, 2013; Phan 
et ai, 2014] and the laboratory [Yamada et ai, 2014]. In the reconnection exhaust, where 
most magnetic energy is released, ion heating takes the form of interpenetrating beams 
[Cowley, 1982; Krauss-Varban and Omidi, 1995; Nakahayashi and Machida, 1997; Hoshino 
et ai, 1998; Gosling et ai, 2005; Lottermoser et ai, 1998; Stark et ai, 2005; Wygant et ai, 
2005; Phan et ai, 2007], which are generated through Fermi rehection in the outhowing, 
contracting magnetic helds. The predicted counterstreaming velocity is twice the exhaust 
velocity CAup in the case of antiparallel reconnection even in the presence of Hall magnetic 
and electric helds [Drake et ai, 2009]. The expected ion temperature increase based 
on such a simple picture is AT* = 0.33 = 0.33 hF. However, in solar wind and 

magnetopause observations the ion temperature increments are signihcantly lower than 
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expected, ATj ~ 0.13 but exhibit the expected scaling with parameters [Drake et ai, 
2009; Phan et ai, 2014], 

The scaling of electron heating is much more challenging to understand because the 
single-pass Fermi reflection yields only a small increase in the electron temperature. Nev¬ 
ertheless, magnetopause observations for electrons yield a similar scaling ATg ~ 0.017 IF 
although with signihcantly less heating compared to the ions [Phan et ai, 2013]. Simula¬ 
tions also yield this scaling [Shay et ai, 2014] although the electron heating mechanism 
remains under debate [Ptaggerty et ai, 2014; Egedal et ai, 2015]. 

Thus, it is important not only to establish the explicit mechanisms for electron and ion 
heating during reconnection but also to determine whether the partition of energy between 
the two species is a universal relation or varies with parameters. We demonstrate here 
through a set of comprehensive computer simulations and analytic methods that the large- 
scale parallel potential that develops within the reconnection exhaust controls and links 
together both electron and ion heating and regulates the partition of released magnetic 
energy. The development of this potential within the exhaust to prevent the escape of hot 
electrons has been well-established [Egedal et ai, 2008] and enables electrons to undergo 
repeated Fermi reflections within the reconnection exhaust. In the present paper we 
identify the mechanism that ultimately limits electron energy gain. The spatial variation 
of the potential propagates outward from the exhaust as a component of a slow shock [Liu 
et ai, 2012]. The electron temperature and the associated shock velocity increase until the 
velocity matches that of the Alfvenic exhaust. At this point electron energy gain through 
Fermi reflection ends since the bounce length of electrons trapped in the exhaust no longer 
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decreases with time as they propagate downstream. At the same time that the potential 
serves to facilitate electron energy gain, it suppresses ion heating: the parallel streaming 
velocity of ions injected into the exhaust from upstream is reduced below the Alfven speed 
by the potential so that the counterstreaming velocity of ions is less than 2caup- Thus, 
the strength of the potential regulates the relative heating of electrons and ions. We show 
that the potential increases with increasing upstream electron temperature Teup and that 
ATg can actually exceed ATj - the partition of electron and ion heating measured in the 
magnetosphere [Eastwood et ai, 2013; Phan et ai, 2013, 2014] and laboratory experiments 
[Yamada et ai, 2014] is not universal. However, the total heating is unaffected by the 
potential and the fraction of magnetic energy converted into thermal energy is constant 
for the simulations performed, with A(Tj + Tg) ~ 0.15 = 0.15 W. Remarkably, 

despite the numerous differences between the simulations and observations, this slope is 
the same as the Phan et al. [2013, 2014] measurement of total heating at the Earth’s 
magnetopause. 

2. Simulations: 

We use the PIC code P3D[Zei/er et ai, 2002] to perform simulations in 2.5 dimensions of 
collisionless antiparallel (no guide held) reconnection. Magnetic held strengths and parti¬ 
cle number densities are normalized to Bq and no, respectively. Lengths are normalized to 
the ion inertial length dio = c/upio at the reference density no, time to the ion cyclotron 
time = {eBo/mic)~^, and velocities to the Alfven speed cao = a/ BQ/{ATTmi no). 
Electric helds and temperatures are normalized to Eq = caqBq/c and Tq = niicj^Q, respec¬ 
tively. In the simulation coordinate system the reconnection outhows are along x and the 
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inflows are along y. Simnlations are performed in a periodic domain with a system size 
of Lx X Ly = 204.8 djo x 102.4 djo, and 100 particles per grid in the inflow region. Simu¬ 
lation parameters, which are given in the table in the Supplementary Material, included 
ion-to-electron mass-ratios of 25 and 100 and a variety of upstream initial temperatures 
and magnetic helds. The initial conditions are a double current sheet [S'da?/ et ai, 2007]. 

A small magnetic perturbation is used to initiate reconnection. Each simulation is 
evolved until reconnection reaches a steady state, and then for analysis purposes during 
this steady period the simulation data is time averaged over 100 particle time steps, which 
is typically on the order of 50 electron plasma wave periods 

3. Overview of electron and ion heating: 

We hrst present an overview of electron and ion heating as measured in the simula¬ 
tions. The temperature of electrons and ions each increase with the distance downstream 
of the X-line in the exhaust until it approaches a constant. This behavior has already 
been discussed in detail for electrons [Shay et ai, 2014] and is discussed more fully in the 
Supplementary Material for the ions. To determine ATj and ATe in a given simulation 
we average T* and Te over a region downstream and then subtract the inflow tempera¬ 
ture. Details of how this average is computed for ions are found in the Supplementary 
Material. In Fig. 1 we present an overview of (a) electron, (b) ion and (c) the total 
temperature increments versus The red triangles correspond to high upstream 

electron temperature Tg/Tj = 9. As expected, the sum of the electron and ion heating 
increments scale with the available magnetic energy per particle with an approximate 
slope of A(Tj -|- Te) ~ 0.15 = 0.15 hE This slope is the same as measured in 
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observations of electron and ion heating in reconnection exhausts at the Earth’s magne¬ 
topause [Phan et ai, 2013, 2014], Surprisingly, however, the individual electron and ion 
temperature increments in Fig. 1 have a larger spread related to the upstream electron 
temperature. The electron heating is generally signihcantly below that of the ions, as in 
the observational data [Phan et ai, 2013]. The exceptions are the runs with high electron 
temperature upstream, which produce enhanced electron heating and reduced ion heating 
with the electron heating signihcantly greater than the ion heating. These simulations 
therefore demonstrate the parameter-dependence of energy partition between electrons 
and ions. In the remainder of the manuscript we explore the mechanisms that control the 
heating of both species, starting with the ions. 

Shown in Figs. 2a-c are the ion parallel temperature Tj||, perpendicular temperature Tij_ 
and total temperature (Tj = |]Tj|| -|- 2Tj_i_] /3). Downstream of the x-line Tj|| increases and 
broadens in the inhow direction to hll the exhaust. The band of Ti± at the midplane of the 
exhaust is produced by the Speiser orbits of the ions [Speiser, 1965; Drake et ai, 2009]. As 
with the electrons, the total ion temperature asymptotes to a constant downstream [Shay 
et ai, 2014]. The underlying mechanism for ion heating well downstream of the x-line was 
outlined by Drake et al. [2009]. In this downstream region, Ej^ = 0 in the reference frame 
moving with the reconnected magnetic held lines. This includes both the reconnecting 
and Hall electric helds as shown in Fig. 3 of Drake et al. [2009]. Within the ion diffusion 
region, however, the strong normal electric held cannot be transformed away [Wygant 
et ai, 2005]. In this moving frame the cold ion population enters the reconnection exhaust 
with a parallel velocity equal to the held line velocity vq. The ions reach the midplane. 
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undergo an energy-conserving reflection, and then travel back out along the held line. The 
reflected population mixes with cold incoming ions creating counter-streaming beams and 
a temperature increment of ATj ^ ATj||/3 ~ rrijUQ/S. In order to test this prediction 
we directly measure the held line velocity Vq ~ —cEz/By, which asymptotes to the ion 
outflow velocity Vix in the downstream region. The prediction of ATj = mjUg/3 is tested 
in Fig. 3a. The points roughly scale with rriiVQ/S, but there are two significant differences 
relative to the theoretical value: (1) There are outlier points leading to a large spread of 
the data, and (2) all of the data points are substantially below the theoretical prediction 
(line of slope = 1 as indicated by the dashed black line). In Fig. 3b, examination of the 
ion distribution function integrated along Vz around (A, Y) = (190, 25.6) reveals that the 
beaming velocities are significantly less than vq. The magnetic held points along y, and 
two field-aligned counterstreaming populations straddle Vy = 0 but well within the region 

\Vy/V0\ < 1 . 

We now show that the reduction in ion heating is a consequence of the large-scale 
potential that confines the hot electrons (see cuts of Tey and n in Fig. 2e) in the exhaust. 
In order to maintain electron force balance along the magnetic held, a large-scale, although 
relativity small magnitude, parallel electric held arises (Fig. 2d). The E\\ Alls the exhaust 
and points away from the midplane. This electric held and associated potential slows 
down inflowing ions leading to a reduced ion beam velocity and a reduced AT*. 

Note that in Fig. 2d there is an inverted E\\ structure straddling the midplane that is 
not to be confused with the larger scale parallel held discussed above. This smaller scale 
parallel electric held is connected with the outer electron diffusion region associated with 
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the super-Alfvenic electron jet [Shay et ai, 2007], and does not couple to the ions, which 
are unmagnetized at these small scales. For that reason, the effect of this electron scale 
parallel electric held is not included in our analysis. 

To calculate the impact of the large scale parallel electric held and the associated poten¬ 
tial on the ions, it is necessary to nnderstand both its amplitude and space-time structure. 
The spatial variation of the potential propagates as a component of the exhaust bound¬ 
ary moving outward away from the midplane. This exhaust boundary takes the form 
of super-slow to sub-slow transition rather than a switch-oh shock because of the strong 
temperature anisotropy that develops in collisionless reconnection[Lm et ai, 2012]. We 
determine this velocity directly from the simulation by calculating the potential 0 by in¬ 
tegrating along the magnetic held. We take 0 = 0 at the X valne of the middle of the 
island {X = 256.4), where the distance along the held line / is also taken to be zero. In 
Fig. 4a we plot 0 versus I and the X-intercept of the held line with the midplane of the 
exhanst, denoted as Xint- Only the portion above the exhaust midplane is shown so that 
the expansion of the white zone with distance downstream measnres the rate of shortening 
of the held line (using the time axis which is dehned by At = AX/vq). The bonndary of 
the white zone parallels the solid line in the white zone, which marks the exhaust velocity 
Vo, so held line shortening is at the velocity uq as expected. The more important resnlt 
of Fig. 4 is that the contours of 0 parallel the boundary of the white zone which means 
that the expansion velocity of the potential is UQ) the same as the shortening rate of the 
held lines. This is a crucial result that will enable us to explicitly calculate ion heating 
and impose limits on electron heating. 
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We analytically calculate the magnitude of 0 from the parallel electric held, which 
follows from electron force balance: 

eE|| = -VfeTeii - Tell Vfolnn + {T^\\ - Te± )VblnT, (1) 

where V;, = (B/T) ■ V. The potential 0, is then given by 0 = — / Ty dl,. Integrating 
Eq. 1 and multiply both side by —1 yields 0 = 0re + (t>n + (t>B) where the subscript 
represents the quantity acted on by the gradient, i.e., 

0n = y(Te||/e) (Vfe \nn)dl. (2) 

In Fig. 4b, these potentials are plotted along the solid black held line shown in Fig. 2d. 
4>Te,4>n, and 0B have diherent constants added to aid in their comparison with 0. 0 
increases from the inhow region to the exhaust, reaching its maximum value just outside 
the midplane. We have found through test particle simulations that as the ions enter the 
exhaust only 0„ signihcantly modihes the ion beam velocity (and therefore ATj). 0^ is 
small. (pTe is signihcant in a narrow region at the edge of the exhaust but because it is so 
localized and because there is a large transverse electric held in this region, test particle 
trajectories provided in the supplementary material reveal that the ions cross this region 
transverse to B and don’t respond to (pTe- The dip in 0 at the midplane of the exhaust 
is similarly unimportant since it only ahects the ion temperature within a narrow region 
that occupies a decreasingly small fraction of the exhaust with distance downstream. 

Thus, <pn has the greatest impact on AT*. To calculate ATj we therefore need to evaluate 
the jump in 0„ across the exhaust. Since Tey is nearly constant across the exhaust (Fig. 2e) 
we can replace it by its average value Tey^ in the integral in Eq. 2. The density varies 
from a minimum rimin at the exhaust boundary to a maximum Ud in the middle of the 
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exhaust so the jump in across the exhaust is given by 


eA0n ^ Te\\d ln{nd/n mm) • 


( 3 ) 


The jump A0„ is marked in the simulation data in Fig. 4b. In Fig. 4c we plot the 
value of A(j)n measured from the simulation against the values from Eq. 3. The agreement 
is excellent. Note the large value of the potential for the simulations with high value of 
upstream T^/Ti (red triangles). We can now extend the model of Drake et al. [2009] to 
include the effect of to obtain a more accurate ion heating prediction. In a frame 
moving with the held line the potential is also unchanging since its outhow velocity is also 
Vo- In this frame the incoming population will be slowed down from the held line velocity 
(uo) to the exhaust beam velocity (vd)- These slower ions mix with incoming ions from 
the other side of the midplane, leading to counterstreaming beams and a temperature 
increment of ATj = miv‘^/3. In the frame of the potential the ion energy is conserved so 
we can calculate Vd directly from ^rriivl — eA^^ = Solving for Vd and substituting 

in ATj, we hnd 


AT,- = 


rrii Vq 


1 - 


2 e A0„ 
rrii vl 


( 4 ) 


In Fig. 4d we insert the measured A^^ from the simulation in Eq. 4 and compare the 
prediction with the measured ion heating in the simulations. The spread in the data is 
markedly reduced compared with that in Fig. 3a: all of the points now straddle a line 
with a slope of 1. Most revealing is the change in position of the Tg/T* = 9 simulations 
which are denoted by red triangles in Figs. 3 and 4. These simulations have large A0„ 
which signihcantly reduces the ion beam velocity and the corresponding ion temperature 
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increment. Thus, electrons, through the self-generated potential have a strong impact on 
ion heating. 

A question remains as to why previous observational studies measure an ion increment 
ATj oc niivl [Drake et ai, 2009; Phan et ai, 2014], even though such a scaling is not 
implied by Eq. 4 due to the presence of the potential A0„. A(f)n depends on both ATe\\, 
and Teup in Eq. 2; the logarithm of the density compression ratio is not expected to 
vary signihcantly with upstream conditions for symmetric reconnection. For a significant 
variation of upstream properties, AT^w has been shown to scale with nii ~ ^ic\up [Phan 
et ai, 2013; Shay et ai, 2014]. Any deviation from the uiivl scaling, therefore, is linked to 
Teup/{,rnic\up) ~ /9enp/2. As long as the electron heating is sufficiently strong compared to 
the upstream temperature, AT^ should dominate, and we recover both the observational 
scaling and the scaling of the black triangles in Fig. 3a. 

Since weak parallel electric fields are impossible to directly measure with in situ satellite 
measurements, the analytic expression for eA(j)n in Eq. 3 can be used to evaluate AT* in 
Eq. 4 to compare with observations. In addition, the prediction can be further simplified 
by using the approximation vq ~ CAup- 

We now discuss the impact of the potential on electron heating. It has been shown that 
the dominant driver of electron heating during anti-parallel reconnection is Fermi reflection 
[Dahlin et al., 2014]. In the absence of scattering, electron energy gain is mostly along 
the local magnetic held. On the other hand, a single Fermi rehection of electrons in the 
reconnection exhaust is not sufficient to drive signihcant electron energy gain. Electrons 
can gain energy through multiple Fermi rehections during multi x-line reconnection [Drake 
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et ai, 2006; Oka et ai, 2010; Drake et ai, 2013] or in a single x-line reconnection as a result 
of the potential 0, which acts to conhne electrons within the reconnection exhaust [Egedal 
et ai, 2008]. What limits the electron temperature within the exhaust and therefore 
the potential (Eq. 2) has not been established. Electrons can continue to gain energy in 
a single exhaust by repeatedly reflecting off of the potential to return to the exhaust core 
for additional Fermi reflections. This behavior is shown by the test particle trajectory in 
Fig 3c, and it is shown in the supplemental material that the reflection is due primarily to 
the potential and not to mirroring. However, electrons lose energy in their reflection from 
the potential (in the frame of the x-line) since the potential is moving outward along the 
magnetic held. The energy gain from Fermi rehection continues to exceed the loss from 
rehection from the potential as long as the expansion velocity is less than vq, the held line 
velocity. Thus, Fig. 4a, which demonstrates that the expansion velocity and held line 
velocity converge downstream, establishes how electron energy gain is limited. The shock 
bounding the reconnection exhaust, as discussed by Liu et al. [2012], carries the potential 
outward along B. The electron temperature increases, increasing the shock velocity, until 
the shock speed reaches Vq and electron heating saturates. 

Here we do not present a complete model of the electron heating during reconnection, 
which requires a full understanding of the dependence of the shock velocity on electron 
and ion temperatures upstream and downstream. Instead we simply note that in the 
limit of low upstream pressure (high upstream mach number) the propagation speed 
of a simple parallel propagating slow shock with a jump in the parallel temperature is 
A{Te + Ti)/mi. Equating this speed to vq = CAup, we hnd A(Te -|-Tj) = 0.25 mi = 
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0.25 W, which is within a factor of two of the simulation and observational hndings of 
0.15 ff^ [Phan et ai, 2013, 2014]. There is signihcant uncertainty in the 0.25 coefficient, 
however, due to the simplistic nature of the shock analysis used to derive it. Nevertheless, 
the basic idea that the electron and ion temperature increments are linked through their 
control of the propagation speed of the shock and associated potential is consistent with 
the results of Fig. Ic. 

4. Conclusions: 

We present the results of PIC simulations of reconnection-driven electron and ion heat¬ 
ing that suggest that the partition of energy gain of the two species is controlled by the 
large-scale potential that develops to prevent hot electrons in the reconnection exhaust 
from escaping along open magnetic field lines. We first show that the relative heating of 
electrons and ions is controlled by the relative magnitudes of the upstream temperatures 
of each species - high upstream electron temperature yields much higher electron than 
ion heating demonstrating that the typical partition of energy seen in space and the labo¬ 
ratory are not universal. We then carry out a detailed study of ion heating and show that 
the potential slows ions injected into the exhaust to values below the Alvenic exhaust flow 
speed. Ion heating therefore can fall well below the characteristic value ATj = rrijUg/S 
predicted by simple Fermi reflection. The scaling of AT^ in the simulations is consistent 
with this theory. The suppression of ion heating becomes very significant for high up¬ 
stream electron temperature when the potential becomes very large. The mechanism by 
which the potential controls electron heating is also discussed. The potential confines 
electrons within the exhaust and enables them to undergo multiple Fermi reflections. The 
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outward propagation of the spatial variation of the conhning potential, which is linked to 
the slow shock that bounds the exhaust, ultimately halts electron energy gain when its 
velocity reaches the exhaust velocity - energy gain through Fermi reflection then balances 
energy loss through reflection off the outward propagating potential. Thus, the electron 
temperature rises until the shock/potential velocity matches the exhaust velocity. The 
potential is therefore the key ingredient that controls both electron and ion heating and 
their relative energy gain. 

An intriguing result is that the total plasma heating {ATfot = AT^ + ATj) in the sim¬ 
ulations is constant with ATtot ~ -15 W, which is consistent with recent magnetospheric 
observations [Phan et ai, 2013, 2014]. Although this is an exciting result, our simulations 
explore only the small parameter regime of symmetric and anti-parallel reconnection. De¬ 
termination of the generality of the ATfot scaling will require a more systematic scaling 
study. Regarding the comparison with satellite observations: On the one hand the fact 
that the observations are of asymmetric reconnection and the simulations are symmetric 
requires some caution during comparison; Clearly, the simulation scaling study should be 
extended to asymmetric reconnection. On the other hand, the fact that asymmetric ob¬ 
servations have such good agreement with symmetric simulations implies that the scaling 
may be a general result, applicable to a wide range of reconnecting systems. 
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Figure 1. Overview of electron and ion heating: (a) ATg, (b) ATj and (c) A(Te + Tj) 
versus mic\^p. The three red triangles have upstream Tg/Ti = 9. The change in total 
temperature appears insensitive to Tgup with A(Te + Ti) ^ 0.15 mic\^p 
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Figure 2. Exhaust structure for a typical simulation with B = 1.0, n = 0.2, Ti = = 

0.25, and mj/rue = 100 (a) Tj||, (b) Ti± and (c) Ti using the same color scale; (d) Spatially 
smoothed Ey with single held line in black. Note that the large scale parallel electric held 
that hlls the exhaust, as opposed to the small scale electric held at the midplane, plays a 
key role in modifying the heating, (e) n, T*, and Tf,\\ along Y at X = 194.5. 
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Figure 3. Test of the basic counterstreaming ion model, (a) ATj versus the theoretical 
prediction m, Ug/3 (the dashed line has a slope of 1). The three red triangles have upstream 
Tg/Tj = 9. (b) Ion distribution function around {X, Y) = (25.6,190) from the simulation 
in Fig. 2. Velocities are normalized to vq (the asymptotic held line velocity), with doted 
red lines showing Vx,y = 0 and dashed yellow lines showing Vy = vq. (c) 2D trajectory of 
a test particle (electron) entering the reconnection exhaust plotted over Bz. The particle 
was initialized upstream with the local ExB velocity, and shows the typical trajectory of 
an electron in the reconnection exhaust. The particle is evolved in the helds from the 
time averaged simulation using the Boris algorithm. 
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Figure 4. (a) The shortening of held lines in the expanding exhaust and the held- 

aligned propagation of the spatial variation of 0 for the same simulation as Fig 2. Shown 

is 0 as a function of distance I along a held line (with / = 0 at the X value of the middle 

of the island (X = 256.9)) and Xint the intercept location of the held line at the midplane 

of the exhaust. The time axis is dehned by At = AXint/vo where Vq is the asymptotic 

held line velocity. 0 is taken to be zero at / = 0. (b) 0, 0Te, 0n, and 0^ versus I along 

the solid black magnetic held line in Fig. 2d with I = 0 dehned as in (a), (c) eA0„ versus 
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Te\\Jn{nd/Umin) ■ (d)ATj versus the predicted temperature increment including the ehect 
of A0„ (Eq. 4). 





















